WHEN IS 0.999 . . . EQUAL TO 1? 
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Abstract. A doubly infinite sum, numerically evaluated at between 0.999 
and 1.001, turns out to have a nice value. 



1. Introduction 

The three dots in the title do not refer to an infinite sequence of 9's, but to digits 
that are increasingly hard to compute. The question is philosophical: how many 
9's do we need to see before we start to believe that the constant we are computing 
is probably equal to 1? In our case, the constant was given by the infinite sum: 

q ~ V H A H k+i - !) m 

where Hj := are the harmonic numbers. 

This question in the title is beyond the scope of a mathematics journal. There 
are, however, mathematical papers proving identities that were discovered because 
numerical computation pointed to a simple answer. The anecdotal evidence then 
accumulates in a misleading manner: when the conjectured identity is false we are 
less likely ever to know. One purpose of the present note is to document a case 
when we were able to evaluate the constant and it turned out not to equal the 
simple guess (in this case, 1). In fact we will prove: 



Theorem 1: 



E 



Hj(Hk+i — 1) 



^ jk{k + l)(j + k) 

-4C(2)-2C(3) + 4C(2)C(3) + 2C(5) (2) 
0.999222... 



where £(n) — 5^ i=1 i ™ denotes the Riemann zeta function. 



The second purpose of this note is to demonstrate a piece of software that, with 
a little human intervention, can find (and prove) this sort of identity. 
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2. Background 



The sum £Q| arose in a paper giving bounds on the run time of the simplex 
algorithm on a polytope known as a Klee-Minty cube |BF04| . The Klee-Minty 
cube is an example of great theoretical importance to analyses of simplex algorithm 
run times. The derivation of the expression Q is, however, of no importance here 
because S appears in |BP04j only as an upper bound, demonstrably not sharp, 
for the leading coefficient, c, of the expected run time. The numerology of S is 
therefore unrelated to the physical origins of c. 

On the other hand, the form of the summation (JTJ does lead one at least to hope 
that an exact value might be derived. In principle, any hypergeometric identity, 
for example, that holds for a finite indefinite or definite sum, may be automatically 
proved via the Wilf-Zeilberger method |P WZ96) . and in fact implementations of 
WZ-type software often can handle summand expressions of the complexity of 
Furthermore, in many cases the WZ-machinery will not only prove but also find 
such an identity, if it exists, given only the left hand side. This, in general, does not 
extend to summations over summands involving no extra parameter: if S were to 
equal 1, that fact would not necessarily be automatically detectable. The relatively 
simple form of the summand, however, gave us hope that summation tricks more 
specific to harmonic series might unlock the problem (the identity ^ l/(n(n + l)) = 
1 stands as a beacon of hope) . Indeed, several telescoping and re-summation tricks 
were initially tried, removing all but one infinite summation in various ways. These 
identities were useful in improving our numerical bounds on S. 

The numerical bounds we had on S were not all that good. It should be noted 
that the harmonic numbers are themselves sums, so the expression is really a 
quadruple sum. This makes it perhaps less surprising that our best rigorous bounds 
were no closer that 10~ 3 . To make a long story short, summing in one variable, 
then using exact values for thousands of row and column sums and an integral 
approximation for the remaining terms, led to our best rigorous bounds, namely 



At this point, although the exact value of S was of no use to us, we felt embarrassed 
to publish numerical bounds on a constant that we suspected was equal to 1. The 
authors of |BP04j then consulted the experts in harmonic summation, who con- 
sulted their computers, and came up with Theorem 1. The remainder of this note 
proves Theorem 1. 



If a sum such as Q has a nice value, it does not necessarily follow that a 
truncated version of the sum has a nice formula. But, if it has a nice closed form, 
automatic methods might succeed in finding it by first computing a recurrence and 
afterwards searching for solutions of the recurrence in which the closed form can 
be represented. Thus we consider a truncated version of the sum, namely, 



0.999197 < S < 1.00093. 



3. Solving the Problem 




b 



a 



that is, the upper limits, instead of being infinite, are taken to be integer variables, 
a and b. In this case, our optimism is rewarded: we will be able to simplify the inner 
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sum so that we can sum a second time if we make some alterations that disappear 
when the upper limit goes to infinity. 

The inner sum. In a first step we compute a closed form evaluation of the inner 
sum 

a 

h(a,k) :=J2f( k >3) 

with f(k,j) := ■ Here we follow the summation principles given in |Sch04| 

that are inspired by IPWZ96| . Note that all the computations are carried out with 
the summation package Sigma in the computer algebra system Mathematica. The 
role of the computer here is to produce equations which we may then rapidly and 
rigorously verify. 

Our first step is to compute for the definite sum h(a,k) the recurrence relation 
k 2 h{a, k)-{k + l)(2fc + 1) h(a, k + 1) 

+ ( » + + 2) „(„, » + 2) = ■*■+*+ *>- <; + ; • (3> 

(k + l)(a + k + l)[a + k + 2) 

with a variation of Zeilberger's creative telescoping trick. 

The way we find this is to guess that there is some r and some constants Co, . . . , c r 
depending on k but not j, and some function g(k,j) which we assume to have a 
relatively simple form, such that a relation 



^2 c sf( k + s > i) = g( k > 3 + 1 )- g( k , i) 



s=0 

holds for all k,j > 1. We ask Sigma to find such a relation for various classes of 
functions g and for r — 1,2,..., until we achieve success with r — 2, Co(fe) = k 2 , 
ci(fc) = -(k + l)(2fc + 1), ca(fc) = (k + l)(k + 2), and g(k,j) = - (fc ffi(£ffii) ■ We 
may then verify the relation 

co(k)f(k, j) + d (fc)/(fc + 1, j) + c 2 (k)f(k + 2, j) = g(k, j + 1) - j) (4) 

by polynomial arithmetic and by using the relation Hj+i = Hj + Summing Q 
over k in the interval {1, . . . , a} proves ©■ 

Next we are in the position to discover and prove that 



kH£-2H k + kH^i +2W {kH a - 1 ^ 1 1^1 
h(a,k) = ——^ — -> :-t> _ , 

v ' 2k 2 k 2 a + i k ^ i ^ a+ i 

i=l i=l j=l J 

(5) 

holds for all a, k > 1; here Hjp = Yli=i F denotes the generalized harmonic 
numbers. 

We do this by asking Sigma to find solutions to J3J for all k > 1 and then to 
plug in the initial conditions, which in this case, are to match the values of h(a, k) 
for k = 1,2; see |Sch04| for further details. Once Sigma has found the expression 
on the right hand side of (JSJ, it is again a finite exercise in polynomial arithmetic 
to verify that it satisfies 0, and that it satisfies the two initial conditions. Since r 
initial conditions uniquely determine the solution, we have proved JHJ. 
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Some terms vanish in the limit. At some point, if we do not winnow out some 
terms that will disappear in the limit, our computer will begin to balk. Luckily it 
is easy to identify some terms that will contribute o(l) to the definite sum as a and 
6 go to infinity, and may therefore be ignored in the evaluation of S. We have, for 
example, the elementary estimates 

1 - ^ . 1 TJ k 2 1 k ^ 2 
lim —7T \^ r = 0, lim — — : = and lim — — ; = 0. 



k 2 ' a + i a^co k ^— ' a + i a— >oo fc -^-^ j ' a + / 

1 i=l i=l j = l •' 



Hence, if we define 



(6) 



6 ^ -Z. ^(fc + 1) ^2 - C'J 

fc=l 

we have 

lim S'(a,b) = S 

a,b — >oo 

by 13 an d ©■ Summarizing, problem |J2J| from above reduces to find and prove 
the identity 

lim S'(a,&) = -4C(2)-2C(3)+4C(2)C(3) + 2C(5). (8) 

a,o — >oo 



The outer sum. Our next step is guided by the fact that we know a few infinite 

sums in which the n th summand is of the form n~ c times a monomial in the har- 

(v) 

monic and generalized harmonic numbers H n and Hn ■ While we cannot solve the 
general problem of summing all such univariate series, it makes sense to attempt to 
manipulate things into this form. Thus we ask Sigma to try to write the summand 
of the right hand side of @ in the form g(a, k) — g(a, k — 1) where the only infinite 
sums that appear in g are of the form described above. The program obligingly 
produces the function g(a, k) = A(a, k) + B(a, k) + C(a, k) where A, B and C are 
given by 

A(a,b) := ^(&H b + AbH b + AHl + 'SbHl + Hl + bHl-QbHW (9) 

2(6 + 1) V 

+2H b H^ + 2bH b HP - 2H (2) - 7bH ( b 2) + H b H ( b 2) + bH b H (2) ^j , 

B(a,b) := -^-Jh^+H^) (10) 
(6+1) v ' 

and 



b TJ b H2 1 b n3 1 b tt rr(2) 

C(a,b) := (^f-DEf -Ef + ^Ef + ^E^" 

i—1 i—1 i—1 i—1 

This time the correctness of the result supplied by Sigma is verified by polynomial 
arithmetic and the definition of Hk without need to appeal to any uniqueness results 
for solutions to recurrences. 
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Summing the relation 

a(k) a(k D H k+1 -lkHl-2H k + kH^+2kHP 
g{k)-g{k 1)- k{k + i) — 

over k in the interval {1, . . . , 6} proves that 

S'{a, b) = A(a, b) + B(a, b) + C(a, b) (12) 
where A, B and C are as in |(?JJ| - (JTTJ. 

Euler sums and multiple £- values. Now we must make good on our supposition 
that the form of g(a, k) leads to sums we can evaluate in terms of the zeta function. 
The limits of A(a, b) and B(a, b) are obvious: 

lim A(a,b)=0 and lim B(a, b) = -4C(2). (13) 

a, fa — too a,b — >oo 

The sums in C(a,b) are handled in two steps, the first being to reduce them to 
multiple (^-values, which can be done for any Euler sum. 

The Euler sum of index (pi < ■ ■ ■ <pk',q), named after |Eul27| (see also |Ber89| ). 
is defined by 

v [/(Pi) fj(P2) rr(Pfe) 

s — V — " ' " (u) 

n=l 

Repeated indices pi = Pi+i are allowed, so the summand may have powers of 
generalized harmonic numbers in the numerator. Each of the summands in C := 
linio^oo C(a, b) is an Euler sum. Define the multiple (- values by 

C(oi, . . . ,a k ) ■ X 



/ , a-i civ. 



«i>— >n k 



A- 



The number k is called the multiplicity of the multiple (Rvalue and a\ + ■ ■ ■ + a k 
is known as the weight. The number a%, . . . ,a k need not be increasing, but it is 
required that a\ > 2, because this is the condition for the sum to be finite. 

One may expand each Hn ** of |Q , an( l m this wa Y the general Euler sum 
becomes the sum of j^ Pl ■ ■ ■ j k Pk n~ q over all (fc + l)-tuples satisfying n > ji, . . . , 
Decomposing according to the set of distinct values among j±, . . . , j k , n produces 
an expression for S pu ,„ lPk - t q that is a linear combination of multiple ^-values with 
weight q + pi + ■ ■ ■ + Pk and multiplicity at most k + 1. In the same way, one sees 
that any product of Euler sums or multiple £- values is itself a linear combination 
of multiple (^-values. 

To illustrate, these ideas, we write C as a sum of products of Euler sums and 
then change these, term by term, into multiple ^-values. From the definitions, 

C = '5'0;2'S'l;2 — Sl;2 — <5l,l;3 + TJ^M,!^ + ^^A 2 ■ 

Of course 5 0;2 = C(2). Next, 

S h 2= £ n- 2 r 1 = E n " 3 + E n- 2 i- X = C(3) + C(2,1). 

n>j>l n>l ti>J>1 
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The remaining values are 

Si,i,3 = 2C(3,1,1)+C(3,2)+2C(4,1) + C(5) 
S ltl ,ifl = 6C(2,1,1,1) + 3C(2,2,1) + 3C(2,1,2) 
+6C(3,1,1) + 3C(3,2) + 3C(4,1)+C(5) 
Si,2; 2 = C(2,2,l) + C(2,l,2) + C(3,2) + C(4,l) + C(5). 
These are derived in the same way as S\-2, for example, 

n>j,k>l 

where the triples decompose into n = j = k,n = j>k,n = k>j,n>j = k,n> 
j > k and n > k > j. Putting this all together gives 

C = (C(2)-1)(C(3)+C(2,1))+C(3,1,1) + C(3,2) + 3C(2,1,1,1) (15) 
+2C(2,2,1) + 2C(2,1,2) + C(2,3). 

The decomposition of products of £- values and multiple (- values into linear combi- 
nations of multiple ^-values is entirely analogous, for instance, 

C(2)C(3) =C(2,3) + C(3,2) + C(5). (16) 

Reducing to ordinary ^-values. Theorem 1 indicates that reduction to single 
^-values is possible. This is worthwhile, because the single ^-values may be more 
easily and accurately obtained than may multiple ^-values. For example, Maple 
computes 500 digits of £(5) in under 10 seconds. 

As it happens, some Euler sums and multiple £- values may be represented as 
sums of products of single (^-values. Others, it appears, may not, though this is 
not proved and no algorithm is known for determining which Euler sums or which 
multiple £- values may be represented in this way. For large k, the number of Q- 
linearly independent multiple ^-values of weight k is conjectured by Zagier to grow 
like 1.32 fe 1FS98I page 17], which outstrips the dimension e '^ of the monomials 
C(fci) • • • ((k m ) in single £-values of total degree k = fci + • • • + k m . However, 
k = 3,4,5,6,7 and 9, the dimensions are known to be equal, and hence every 
multiple £-value of these weights is a polynomial over Q in single £- values. 

Among the known relations, two classes suffice for our problem, namely the sum 
and duality relations. The sum relation of degree k says that the sum of all multiple 
^-values of fixed multiplicity m and weight k is equal to C(&)- F° r k — 3,5 and 
m = 2, these are 

C(3) = C(2,l) 

C(5) = C(4,l) + C(3,2) + C(2,3). 

The sum relations were conjectured by Moen |Hof92l Section 3], known in special 
cases to Euler Eul27j , and proved by |Gra97j . 

Secondly, one has the duality theorem, conjectured by Ho f92| and proved by 
|Zag94| . These imply equality between C(Pi> ■ • ■ :Pk) and ((qi, ...,©) of equal 
weights when the index sequences are related by an operation akin to transpos- 
ing a partition. The exact definition is given in, e.g., |FS98I page 29] but all we 
will need to here is to represent the multiple ^-values of weight 5 and multiplicity 3 
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and 4 in terms of those of multiplicity 2: 

C(2,2,l) = C(3,2) 

C(2, 1,2) = C(2,3) 

C(3,l,l) = C(4,l) 

C(2, 1,1,1) = C(5). 

Using the duality relations to eliminate all terms of multiplicity greater than 2 
from ljl5|) . and then using (|16J) . gives 

C = (C(2) - 1)(C(3) + C(2, 1)) + C(4, 1) + 3C(2)C(3) . 

The sum relations now yield 

C = 4C(2)C(3)-2C(3) + 2C(5), 

which finishes the proof of Theorem 1. 

We remark that one more relation is needed to reduce all multiple ^-values of 
weight 5 to polynomials in single ^-values, namely 

C(3,2)=3C(2)C(3)-y C (5) J (17) 

but that our expression for C can be simplified in terms of single £-values with- 
out Q17[l. This identity 117|l can be produced by a formula in BBG9 5 which can 
transform any double C-value with odd weight to a polynomial in single C-values. 
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